library(tidyverse)
library(tuneR)
library(devtools)
library(ggplot2)
library(tsfeatures)
library(class)
library(cvTools)
library(randomForest)
library(GEOquery) 
library(R.utils)
library(reshape2)
library(limma)
library(dplyr)
library(e1071)
library(DT)
library(viridis)
library(plotly)
library(scales)
gse36059 = getGEO("GSE36059")[[1]]
Found 1 file(s)
GSE36059_series_matrix.txt.gz
Using locally cached version: C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GSE36059_series_matrix.txt.gz
Warning in for (getMethodName in getMethodNames) { :
  closing unused connection 5 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE36nnn/GSE36059/matrix/)
Warning in .Internal(lapply(X, FUN)) :
  closing unused connection 4 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE36nnn/GSE36059/matrix/)
Warning in .Internal(lapply(X, FUN)) :
  closing unused connection 3 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE48nnn/GSE48581/matrix/)
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse48581 = getGEO("GSE48581")[[1]]
Found 1 file(s)
GSE48581_series_matrix.txt.gz
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse129166 = getGEO("GSE129166")[[1]]
Found 1 file(s)
GSE129166_series_matrix.txt.gz
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse36059_f = fData(gse36059)
gse36059_f
gse48581_f = fData(gse48581)
gse48581_f
gse129166_f = fData(gse129166)
gse129166_f
gse36059_p = pData(gse36059)
gse36059_p
gse48581_p = pData(gse48581)
gse48581_p
gse129166_p = pData(gse129166)
gse129166_p
#Remove Nephrectomy outcomes
gse36059_p = gse36059_p[!(gse36059_p$characteristics_ch1=="diagnosis: Nephrectomy"),]

gse48581_p = gse48581_p[!(gse48581_p$characteristics_ch1=="diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): nephrectomy"),]
#Encodes stable as 0 and rejecting as 1
gse36059_p$diagnosis = ifelse(gse36059_p$characteristics_ch1 == "diagnosis: non-rejecting", 0, 1)
gse48581_p$diagnosis = ifelse(gse48581_p$characteristics_ch1.1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): non-rejecting", 0, 1)
gse129166_p$diagnosis = ifelse((gse129166_p$characteristics_ch1.1 == "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 0, 1)
#Encodes stable as 0, ABMR as 1, TCMR as 2, Mixed as 3
gse36059_p$exact_diagnosis = ifelse(gse36059_p$characteristics_ch1 == "diagnosis: non-rejecting", 0, ifelse(gse36059_p$characteristics_ch1 =="diagnosis: ABMR", 1, ifelse(gse36059_p$characteristics_ch1 == "diagnosis: TCMR", 2, 3)))

gse48581_p$exact_diagnosis = ifelse(gse48581_p$characteristics_ch1.1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): non-rejecting", 0, ifelse(gse48581_p$characteristics_ch1.1 =="diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): ABMR", 1, ifelse(gse48581_p$characteristics_ch1.1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): TCMR", 2, 3)))

gse129166_p$exact_diagnosis = ifelse((gse129166_p$characteristics_ch1.1 == "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 0, ifelse((gse129166_p$characteristics_ch1.1 != "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 != "abmr (no: 0_Yes:1): 0"), 3, ifelse((gse129166_p$characteristics_ch1.1 != "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 2, 1)))
gse36059_p
gse48581_p
gse129166_p

Main outcomes are non-rejecting, TCMR (acute T-cell–mediated rejection), ABMR (anti-donor antibody-mediated rejection), and MIXED. Also has nephrectomy.

gse36059_ex = data.frame(t(exprs(gse36059)))
gse36059_ex
gse48581_ex = data.frame(t(exprs(gse48581)))
gse48581_ex
gse129166_ex = data.frame(t(exprs(gse129166)))
gse129166_ex
start = 1
stop = 100
boxplot(gse36059_ex[start:stop])

boxplot(gse48581_ex[start:stop])

boxplot(gse129166_ex[start:stop])

LS0tDQp0aXRsZTogIkxpYW0tUmVzZWFyY2gtMi1JbnRlZ3JhdGlvbiINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCmBgYHtyfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KHR1bmVSKQ0KbGlicmFyeShkZXZ0b29scykNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkodHNmZWF0dXJlcykNCmxpYnJhcnkoY2xhc3MpDQpsaWJyYXJ5KGN2VG9vbHMpDQpsaWJyYXJ5KHJhbmRvbUZvcmVzdCkNCmxpYnJhcnkoR0VPcXVlcnkpIA0KbGlicmFyeShSLnV0aWxzKQ0KbGlicmFyeShyZXNoYXBlMikNCmxpYnJhcnkobGltbWEpDQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShlMTA3MSkNCmxpYnJhcnkoRFQpDQpsaWJyYXJ5KHZpcmlkaXMpDQpsaWJyYXJ5KHBsb3RseSkNCmxpYnJhcnkoc2NhbGVzKQ0KYGBgDQoNCmBgYHtyfQ0KZ3NlMzYwNTkgPSBnZXRHRU8oIkdTRTM2MDU5IilbWzFdXQ0KZ3NlNDg1ODEgPSBnZXRHRU8oIkdTRTQ4NTgxIilbWzFdXQ0KZ3NlMTI5MTY2ID0gZ2V0R0VPKCJHU0UxMjkxNjYiKVtbMV1dDQpgYGANCg0KDQoNCmBgYHtyfQ0KZ3NlMzYwNTlfZiA9IGZEYXRhKGdzZTM2MDU5KQ0KZ3NlMzYwNTlfZg0KZ3NlNDg1ODFfZiA9IGZEYXRhKGdzZTQ4NTgxKQ0KZ3NlNDg1ODFfZg0KZ3NlMTI5MTY2X2YgPSBmRGF0YShnc2UxMjkxNjYpDQpnc2UxMjkxNjZfZg0KYGBgDQoNCmBgYHtyfQ0KZ3NlMzYwNTlfcCA9IHBEYXRhKGdzZTM2MDU5KQ0KZ3NlMzYwNTlfcA0KZ3NlNDg1ODFfcCA9IHBEYXRhKGdzZTQ4NTgxKQ0KZ3NlNDg1ODFfcA0KZ3NlMTI5MTY2X3AgPSBwRGF0YShnc2UxMjkxNjYpDQpnc2UxMjkxNjZfcA0KYGBgDQpgYGB7cn0NCiNSZW1vdmUgTmVwaHJlY3RvbXkgb3V0Y29tZXMNCmdzZTM2MDU5X3AgPSBnc2UzNjA1OV9wWyEoZ3NlMzYwNTlfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxPT0iZGlhZ25vc2lzOiBOZXBocmVjdG9teSIpLF0NCg0KZ3NlNDg1ODFfcCA9IGdzZTQ4NTgxX3BbIShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDE9PSJkaWFnbm9zaXMgKHRjbXIsIGFibXIsIG1peGVkLCBub24tcmVqZWN0aW5nLCBuZXBocmVjdG9teSk6IG5lcGhyZWN0b215IiksXQ0KYGBgDQoNCg0KYGBge3J9DQojRW5jb2RlcyBzdGFibGUgYXMgMCBhbmQgcmVqZWN0aW5nIGFzIDENCmdzZTM2MDU5X3AkZGlhZ25vc2lzID0gaWZlbHNlKGdzZTM2MDU5X3AkY2hhcmFjdGVyaXN0aWNzX2NoMSA9PSAiZGlhZ25vc2lzOiBub24tcmVqZWN0aW5nIiwgMCwgMSkNCmdzZTQ4NTgxX3AkZGlhZ25vc2lzID0gaWZlbHNlKGdzZTQ4NTgxX3AkY2hhcmFjdGVyaXN0aWNzX2NoMS4xID09ICJkaWFnbm9zaXMgKHRjbXIsIGFibXIsIG1peGVkLCBub24tcmVqZWN0aW5nLCBuZXBocmVjdG9teSk6IG5vbi1yZWplY3RpbmciLCAwLCAxKQ0KZ3NlMTI5MTY2X3AkZGlhZ25vc2lzID0gaWZlbHNlKChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjEgPT0gInRjbXIgKG5vOiAwX2JvcmRlcmxpbmU6MV9UQ01SOjIpOiAwIikgJiAoZ3NlMTI5MTY2X3AkY2hhcmFjdGVyaXN0aWNzX2NoMS4yID09ICJhYm1yIChubzogMF9ZZXM6MSk6IDAiKSwgMCwgMSkNCmBgYA0KDQoNCmBgYHtyfQ0KI0VuY29kZXMgc3RhYmxlIGFzIDAsIEFCTVIgYXMgMSwgVENNUiBhcyAyLCBNaXhlZCBhcyAzDQpnc2UzNjA1OV9wJGV4YWN0X2RpYWdub3NpcyA9IGlmZWxzZShnc2UzNjA1OV9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0gImRpYWdub3Npczogbm9uLXJlamVjdGluZyIsIDAsIGlmZWxzZShnc2UzNjA1OV9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0iZGlhZ25vc2lzOiBBQk1SIiwgMSwgaWZlbHNlKGdzZTM2MDU5X3AkY2hhcmFjdGVyaXN0aWNzX2NoMSA9PSAiZGlhZ25vc2lzOiBUQ01SIiwgMiwgMykpKQ0KDQpnc2U0ODU4MV9wJGV4YWN0X2RpYWdub3NpcyA9IGlmZWxzZShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSA9PSAiZGlhZ25vc2lzICh0Y21yLCBhYm1yLCBtaXhlZCwgbm9uLXJlamVjdGluZywgbmVwaHJlY3RvbXkpOiBub24tcmVqZWN0aW5nIiwgMCwgaWZlbHNlKGdzZTQ4NTgxX3AkY2hhcmFjdGVyaXN0aWNzX2NoMS4xID09ImRpYWdub3NpcyAodGNtciwgYWJtciwgbWl4ZWQsIG5vbi1yZWplY3RpbmcsIG5lcGhyZWN0b215KTogQUJNUiIsIDEsIGlmZWxzZShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSA9PSAiZGlhZ25vc2lzICh0Y21yLCBhYm1yLCBtaXhlZCwgbm9uLXJlamVjdGluZywgbmVwaHJlY3RvbXkpOiBUQ01SIiwgMiwgMykpKQ0KDQpnc2UxMjkxNjZfcCRleGFjdF9kaWFnbm9zaXMgPSBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSA9PSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgPT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAwLCBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSAhPSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgIT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAzLCBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSAhPSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgPT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAyLCAxKSkpDQoNCmBgYA0KDQoNCmBgYHtyfQ0KZ3NlMzYwNTlfcA0KZ3NlNDg1ODFfcA0KZ3NlMTI5MTY2X3ANCmBgYA0KTWFpbiBvdXRjb21lcyBhcmUgbm9uLXJlamVjdGluZywgVENNUiAoYWN1dGUgVC1jZWxs4oCTbWVkaWF0ZWQgcmVqZWN0aW9uKSwgQUJNUiAoYW50aS1kb25vciBhbnRpYm9keS1tZWRpYXRlZCByZWplY3Rpb24pLCBhbmQgTUlYRUQuIEFsc28gaGFzIG5lcGhyZWN0b215Lg0KDQoNCmBgYHtyfQ0KZ3NlMzYwNTlfZXggPSBkYXRhLmZyYW1lKHQoZXhwcnMoZ3NlMzYwNTkpKSkNCmdzZTM2MDU5X2V4DQpnc2U0ODU4MV9leCA9IGRhdGEuZnJhbWUodChleHBycyhnc2U0ODU4MSkpKQ0KZ3NlNDg1ODFfZXgNCmdzZTEyOTE2Nl9leCA9IGRhdGEuZnJhbWUodChleHBycyhnc2UxMjkxNjYpKSkNCmdzZTEyOTE2Nl9leA0KYGBgDQpgYGB7cn0NCnN0YXJ0ID0gMQ0Kc3RvcCA9IDEwMA0KYm94cGxvdChnc2UzNjA1OV9leFtzdGFydDpzdG9wXSkNCmJveHBsb3QoZ3NlNDg1ODFfZXhbc3RhcnQ6c3RvcF0pDQpib3hwbG90KGdzZTEyOTE2Nl9leFtzdGFydDpzdG9wXSkNCmBgYA0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQo=